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ABSTRACT 

Afterglow light curves are constructed analytically for realistic gamma-ray 
burst remnants decelerating in either a homogeneous interstellar medium or a 
stellar wind environment, taking into account the radiative loss of the blast wave, 
which affects the temporal behaviors significantly. Inverse Compton scattering, 
which plays an important role when the energy equipartition factor e e of electrons 
is much larger than that of the magnetic field (eg), is considered. The inverse 
Compton effect prolongs the fast-cooling phase markedly, during which the rela- 
tivistic shock is semi-radiative and the radiation efficiency is approximately con- 
stant, e = e e . It is further shown that the shock is still semi-radiative for quite a 
long time after it transits into the slow-cooling phase, because of a slow decreas- 
ing rate of the radiation efficiency of electrons. The temporal decaying index of 
the X-ray afterglow light curve in this semi-radiative phase is (3p — 2 + 2e) /(4 — e) 
in the interstellar medium case, and [3p — 2 — (p — 2)e]/2(2 — e) in the stellar 
wind case, where p is the distribution index of the shock-accelerated electrons. 
Taking p = 2.2 — 2.3 as implied from common shock acceleration mechanism, 
and assuming e e ~ 1/3, the temporal index is more consistent with the observed 
(ax) ~ 1-3 than the commonly used adiabatic one. The observability of the 
inverse Compton component in soft X-ray afterglows is also investigated. To 
manifest as a bump or even dominant in the X-ray afterglows during the rela- 
tivistic stage, it is required that the density should be larger than ~ 1 — 10 cm -3 
in the interstellar medium case, or the wind parameter A* should be larger than 
~ 1 in the stellar wind case. 



Subject headings: gamma rays: bursts — hydrodynamics — stars: mass loss - 
ISM: jets and outflows — radiation mechanisms: nonthermal — relativity 
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1. Introduction 

In the past seven years, tremendous progress in understanding the cosmic gamma- 
ray bursts (GRBs) has been achieved from the detections of long-lived GRB embers (or 
afterglows) in low frequency bands (for reviews see van Paradijs, Kouveliotou, & Wijers 
2000; Cheng & Lu 2001; Meszaros 2002; Zhang & Meszaros 2004; Piran 2004). The simple 
standard shock model (Meszaros & Rees 1997; Sari, Piran, & Narayan 1998) responsible for 
the afterglow has been successfully established (Wijers, Meszaros, & Rees 1997; Vietri 1997; 
Waxman 1997; Galama et al. 1998). 

As more and more afterglows being detected, the energetics of GRB remnants and shock 
microphysics have been inferred within the context of the standard model. Although the 
standard model are in rough agreement with afterglow observations, problems still exist 
since the standard scenario is oversimplified at least in two aspects. First, in the standard 
afterglow model the relativistic shock is usually assumed to be quasi-adiabatic. However, 
the shock in fact may be partially radiative. Fittings to observed afterglows reveal that the 
shock imparts an equipartition amount of energy into electrons, which is responsible for both 
the afterglow emission and the shock energy loss (Panaitescu & Kumar 2001, 2002). The 
temporal evolution of the shock energy would affect the estimation of the GRB energetics 
from late time afterglows, as well as the profiles of afterglow light curves (Lloyd-Ronning 
& Zhang 2004; Bottcher & Dermer 2000). Second, multi- wavelength fittings to several 
afterglows also indicate that post-shock energy density imparted to electrons is statistically 
much larger than that of the magnetic fields (Panaitescu & Kumar 2001). It implies that 
the inverse Compton (IC) scattering plays an important role in GRB afterglows. The IC 
scattering has two effects. One is to enhance the cooling rate of the shock-accelerated 
electrons and hence delay the transition from early fast-cooling phase to late slow-cooling 
phase. It finally influences the observed flux density. Another effect is to cause a high energy 
spectral component typically above the soft X-ray band. These IC effects have been taken 
into account by numerous authors (Waxman 1997; Wei & Lu 1998, 2000; Panaitescu & 
Kumar 2000; Dermer, Bottcher, & Chiang 2000; Sari & Esin 2001; Bjornsson 2001; Wang, 
Dai, k Lu 2001; Zhang & Meszaros 2001; Li, Dai, & Lu 2002). 

As for the circum-burst environment, in the standard model it is once assumed that the 
surroundings of GRBs are homogeneous interstellar medium (ISM). Over the past several 
years, a lot of evidence has been collected linking GRBs to the core collapse of massive 
stars (Woosley 1993; Paczyhski 1998). The most important evidence came from the direct 
association between GRB 030329 and the supernova SN 2003dh (Stanek et al. 2003; Hjorth 
et al. 2003), as well as the previous tentative association between GRB 980425 and SN 
1998bw (Kulkarni et al. 1998; Galama et al. 1998). These associated supernovae have been 
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confirmed to be type Ib/c SNe. The progenitors of type Ib/c SNe are commonly recognized 
as massive Wolf-Rayet stars. During their whole life, the massive progenitors eject their 
envelope material into their surroundings through line pressure and thus the stellar wind 
environments are formed. This means that the circum-burst medium for GRB afterglows 
may be the stellar wind (Dai & Lu 1998; Meszaros, Rees, & Wijers 1998; Chevalier & Li 
1999, 2000; Panaitescu & Kumar 2000). 

In this paper, we study the afterglow properties of realistic GRB shocks, considering 
the effect of energy losses. The circum-burst environment is assumed to be either the ISM- 
type or the stellar wind type. We present an analytical solution for the realistic blast wave 
during the fast-cooling phase of GRB afterglows in §2. This semi-radiative hydrodynamics is 
applied to the late slow-cooling phase with quite reasonable argument in §3. Constraints on 
the IC components in the soft X-ray afterglows are given in both sections. In §4 we illustrate 
typical analytical light curves for the realistic model in detail. Conclusions and discussion 
are presented in §5. 



The realistic model for GRB remnants has been extensively investigated in the past 
few years (Huang, Dai, & Lu 1999; Huang et al. 2000). It has been shown that this model 
is correct for both adiabatic and radiative fireballs, and in both ultra-relativistic and non- 
relativistic phases. The basic hydrodynamic equation of this model can be derived as follows. 
In the fixed frame, which is rest to the circum-burst environment, the total kinetic energy of 
the fireball is E K = (7 - l)(M cj + M sw )c 2 + (1 - e)^U, where 7 is the Lorentz factor of the 
blast wave, M e j is the initial mass of the blast wave ejected from the central engine, M sw is 
the mass of the swept-up ambient medium, c is the speed of light, and e is the total radiation 
efficiency (Huang et al. 1999). In the comoving frame of the blast wave, the total internal 
energy instantaneously heated by the shock is C/ = (7 — l)M sw c 2 , which is implied from the 
relativistic jump conditions (Blandford & McKee 1976). The differential loss of the kinetic 
energy Ek, when the blast wave sweeps up an infinitesimal mass dM sw , can be formulated 
as 



Assuming a constant e and inserting the expression for U, it is then easy to obtain the 
hydrodynamic equation of the realistic model (Huang et al. 1999, 2000) 



2. The early fast-cooling phase 



d[(7-l)(M cj + M sw ) C 2 + (l- e ) 7 t/] 



e7 (7-l)dM sw c 2 . 



(1) 



d7 




(2) 



dM 



sw 



M ej +eM sw + 2(l-e) 7 M sw ' 
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Feng et al. (2002) have relaxed the assumption of a constant e, and found that the results 
differ little from the above equation. As we show later, the epoch for a constant radiation 
efficiency will not end just at the transition from the fast cooling phase to the slow cooling 
phase of the fireball evolution. In fact, it would last for a time much longer than that 
transition time, because the radiation efficiency of the shock-accelerated electrons in the 
slow cooling phase decreases very slowly for typical values of the index of electron energy 
distribution, e.g. p ~ 2.2 indicated both from observations of the afterglow spectra, and 
from the shock acceleration theory (Achterberg et al. 2001). 

Throughout this paper, we focus on the early epoch when the afterglow of a relativistic 
jet is spherical-like, which requires the Lorentz factor of the jet (7) being larger than the 

ATT n 

inverse of the half-opening angle. The swept-up mass is given by M sw = -m„nR , where 

3 — k 

m p is the proton mass, and the ambient density is 



where k = with n = A = const for the ISM case, and k — 2 with A = 3 x lO 35 ^ cm -1 
for the stellar wind case. Such a blast wave begins to decelerate at the radius -Rdec, when 
the swept-up mass reaches M e j/7o, where 70 is the initial Lorentz factor. The corresponding 
decelerating time measured by an observer is tdcc = -Rdcc(l + z )/2lo c , here z is the cosmo- 
logical redshift of the GRB. We neglect the effect of reverse shocks in the early afterglow for 
simplicity. 

At early times, electrons cool rapidly. The blast wave is therefore semi-radiative with 
a constant radiation efficiency e = e e . The typical energy equipartition factor of electrons is 
e e ~ 1/3. Equation (2) can then be analytically integrated by neglecting the first two terms 
in the denominator at the right side when t > tdcc- The scaling laws for the hydrodynamics 
are 7 2 oc R~ m and R oc t 1 /( m+1 ), where the hydrodynamic self-similarity index 



and t is the observed time since the burst. The e term in the denominator in equation 
(4) shows the deviation of the hydrodynamics of a semi-radiative blast wave from that of 
an adiabatic one of Blandford & McKee (1976). Since the isotropic-equivalent energy E is 
proportional to M sw 7 2 , it decreases as 



where E& ec is the initial isotropic energy at i?dec- The minimum Lorentz factor of the shock- 
accelerated electrons is evaluated by 



n = AR~ k , 



(3) 



E = £dcc(-V (3 ~ fcM4 - fc - e) , 



(5) 



^dec 




t ™ < \-m/[2(m+l)l 
f-ricc 



dec 



(6) 
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p-2 

where Ci/6 = 6 and m e is the electron mass. Conventionally, assuming a constant 

p — 1 

fraction of e# of the post-shock thermal energy density contained in post-shock magnetic 
fields, the magnetic field intensity in the comoving frame is 

B = B dcc (—r im+k ^ 2{m+ ^, (7) 

^dec 

where the initial value is -Bdcc = (32neBAR^m p c 2 ) 1 ^ 2r yo. The maximum Lorentz factor of 
the shock-accelerated electrons, 7 eimax ~ 10 8 (S/G)~ 1 / 2 , is obtained by assuming that the 
acceleration timescale, which is typically the gyration period in the magnetic field, equals 
the hydrodynamical timescale. The cooling Lorentz factor of electrons is determined by 
considering both synchrotron radiation and inverse Compton scattering, i.e. (Sari et al. 
1998; Panaitescu & Kumar 2000) 

= 6nm e c(l + z) t ,( m+2 k-2)/\2( m +i)} /on 
lc a T 5 2 ec7 otdec(l + mdcc ' K> 

where o"t is the Thomson cross section, and the Compton parameter Y = (— ) 1//2 ^> 1 

is a constant during the fast cooling phase (Sari & Esin 2001). The transition from the 
fast-cooling phase to the slow-cooling phase happens at t cm , when 7 C = 7 eiin i n . Combining 
equations (4) — (8) with the definition of the deceleration time, we obtain 



1/2 

_ (1 + Z)(7 T ^(3-fc)/2 t (3 - fc)^c m i (2^fc)/2/ 2m p ,3/4,l/4 A l/2^4-fc / Q \ 

tcm ~ 2c AaT [ 4nm n ci 1 l 3m P e =1/6/ ■ 



fe.mrn- 



Here the subscript "cm" denotes the physical quantity at the time when 7 C = 7 e 
Throughout this paper we are especially interested in the usual case of e e ^> in which 
the inverse Compton scattering has a dominant effect on the evaluation of 7 C . We denote 
the isotropic energy at t cm as E cm = E(t cm ). Note that t cm is independent of the radiation 
efficiency e and the initial Lorentz factor 7 . The value of t cm can be further calculated to 
be 

t = f 0.30(1 + z)<? _ . 5 eB,-2. 5 g /6 E cm , 53 n day, ISM, 
" 0.58(1 + 2)e^ 05 eJ/ 2 _ 25 Ci/6^* day, wind. 

Here we adopt the conventional definition of Q = Q X 10 X . The radius of the blast wave at 



tcm IS 



_ 2m 



p ^3/4.1/4^1/2 r (3 - K)Ecm 1/2 1/2 nn 
~ e e e B Si/6 L 7— o J °T > U 1 , 



3m e 6 47rmpC 2 



or equivalently 



, 3.98 x 10" f f o.s^i^B'/^ cm, ISM, 
= ' 2.30 x cm, wind. < l2 > 
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The evolution of the radius is R = R cm (——Y 1 k ^- The magnetic field at t cm is 

tr. 



D _ K 3 U/4 m e c4 3m e ^^9/ 8 ^l/8^-3/4 r (3 ~ k)E cm ] _ 1/A 

& cm - o{—) „ e e e B c, 1/6 [ j , ^roj 



'2tt ; e 3 2m p e Sl / 6 1 Anm p c 2 

or numerically 

-9/8 1/8 a-3/4 TTi-1/4 



0.35e e _ 5 e B _ 2 gC^g E cm b3 G, ISM, 

-Dcm — ^ _g/8 1/8 >-3/4 p -l/4 n . , l i4 J 

- 46e e,-0.5 e B,-2. 5 Cl/6 ^cm,53 G > wmd - 



The magnetic field evolves as B = B cm ( J_)-( 3 -* e )/[2(4-fc- e )]. 



2.1. Properties of the synchrotron radiation 

The characteristic synchrotron frequencies corresponding to the j c , 7e im i n and 7 e , m ax 

electrons are denoted as u c , u rn and z/m respectively. They can be easily calculated according 
eB 

to v = — r , with e being the electron charge. The peak flux density of the 

27r(l + z)m e c 

afterglow is F^ max = ? ^^ (l + z), where N e = -AR 3 ~ k is the total number of 

4:7tD l 3 — k 

shocked accelerated electrons, P^max = — 7-B is the peak spectral power of a single 

3e 

electron and is the luminosity distance (Sari et al. 1998). This peak flux density is at the 
cooling frequency v c in the fast-cooling phase, and the flux density at v c would be reduced 
if the synchrotron-self-absorption (SSA) frequency v a is above v c . 

It is convenient to re-scale the physical quantities to the values at the time t cm , because 
physical variables such as 7 C and / ~f m at t cm are independent of e. The minimum electron 
Lorentz factor equals to the cooling Lorentz factor at t cm , 7e )Cm = 7 e ,min(^cm) — 7c(^cm)> 

_ 1 2mp v( fc -l)/2 p (3fc-l)/8 Jfc-3)/8 Ak+\)/A / 4 (3-fc)/2x-i/ 2 r (g ~ fcjjjcm i (fe-l)/4 /-, r\ 

7e ' cm ~4 l 3W £e Cfl Cl / 6 1 T j [ 47rm pC 2 1 ' U5j 
which can be evaluated to be 

f 1 ^3fi^ 1/8 f ~ 3/8 r 1 /4/r- 1 /4_-l/2 TSM 

I iooDe ei _ 05 e Bi _ 25 (, 1/6 £/ cm!53 n , ism, 

7e,cm — S 5/8 -1/8 .3/4 pl/4 ,-1/2 . , \ W ) 

[ 849e e ;_ a5 e B; i 25 C 1 ; 6 £ cm ,53^* > wmd - 
The typical frequency v m also equals to the cooling frequency u c at t cm , i.e. z/ cm = ^ m (t cm ) = 
Vc(t C m), which is 

_ 1 , 3 5/4 m e C 3 / 2rnp .( 3fe _ 7 ) /2 ( 9fe _2Q)/8 (3fc-8)/8^(3fe-4)/4 

^cm - i«V „J ^2/1 , ^U™ ^ e e e B Si/6 



x ( Aat k)/2 ) - 3/2 [ (3 ~ fc)i ! cm ] (3fc " 4)/4 , (17) 



47rm p c 2 
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and can be further deduced to be 

J 3.65 x 10 13 (1 + z)-X,%e- B ^ 5 CylE^ 53 n-^ 2 Hz, ISM, 
" I 8.10 x 10»(1 + zy^^^^' 2 Hz, wiud. (18) 

The maximum frequency of the synchrotron radiation at t cm is 

/, \ _ ( ^-3 X 10 25 (l + z) 1 e e ,- / 0.5 e _B,-2.5Ci/6-^'cm^3'-" " , ,, . 

VM\Lcm) — \ n . ln9R/l , wl 5/8 -1/8 .3/4^1/4 .-1/2 TT . , l iy J 



4.3 x 10 25 (1 + z)-X%etX^,tE2%n- 1 ' 2 Hz, ISM, 

2.4 x 10 25 (1 + ^-^fo.s^LC^s^" 172 Hz, wind, 



which corresponds to ~ 100 GeV photons. It ensures that the synchrotron spectrum can be 
extrapolated to very high energy band, as will be useful in the next subsection. The peak 
flux density of the synchrotron radiation at t cm is 

p (f \ _ ,2?r 3/4 4m e c 2 (l + z) 3m e (k _ 1)/2 3k/8 {4 _. k)/8 k/4 
"' maxlcmj " l 3 j (3 -*;)£>£ l 2m/ e 

x(A4 3 - fc)/2 ) 1 / 2 [ ^— ]( 4 - fc )/ 4 , (20) 
A7rm p c 2 

which is numerically expressed as 



p u i_J 44 (1 + ^'B^^-B^n 1 ' 2 ^ mjy, ISM, 
WU) " < 104(1 + ^Ct^C/W^Xls mJy, wind. U ' ' 



We obtain the temporal evolutions of these characteristic frequencies and the peak flux 
density during the fast-cooling phase as follows, 

_ f_L^(3fe-4)/[2(m+l)] _ f J_^-(4m+fc)/[2(m+l)] 

"C "cm \. J 5 "m "cm \,J 5 

™ w 

^ = ^M(t Cm )(-^)- m/[2(m+1)1 , i^nax = ^, m ax(t cm )(-^) (6 - 3fc - 2m)/[2(m+1)] . (22) 

The synchrotron-self-absorption frequency in the fast cooling phase can be evaluated by 

r c eni? l3/5 

(3 - fc) r2 „x 
r c enit 11/3 .„ v ; 

^[7^ TT^tJ > lf < < I'm, 



■(3-*) B 7| 



jr> _|_ 2 

where c ~ 10.4 — (see the appendix of Wu et al. 2003). When t < t cm , the distribution 

p + 2/3 

of the cooled electron with 7 C < 7 C < 7e,min has p = 2 and the resulting coefficient c = 15.6. 
For 2 < p < 3, the value of Co is nearly a constant, ~ 15, which is about 3 times larger than 
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the coefficient in equation (52) of Panaitescu & Kumar (2000). The SSA frequency can be 
determined straightforwardly by 

v a = min{z/ aj< , v a> } (24) 

without judging whether u a < u c or not. The case for u a > u m can be neglected. The 
numerical expression for z/ a < is 

4.75(1 + .)- 1 e -io.5^%. 5 Cr / 1 6 <53n 3 / 5 (£-)- (10+8e)/[5(4 - £)] GHz, ISM, 
^ ~ 20.9(1 + zYXX^l^T ^S^ 75 ?^)^ 16 - 10 ^^ 2 "^ GHz, wind, (25) 
while the expression for z/ a > is 

2.53 x 10"(1 + ^-^-^ae-^LC/e^g^- 17 ^^)" 27 ^^ ^ EM, 

V = I tcm C26) 

"■ S 2.96 x 10"(1 + ^-^-/o^e^LCvfA- 1 /^-^)- 2 / 3 Hz, wind. 

For the ISM case, the transition from initially u a = v a > to the later v a = v a < happens when 

v a = v c , 

f 6.03 x 10- 51 e^ 1 5 e^_ 2 . 5 ^ i53 n 22 - 75 t cm , if e = 0.1, 
tac I 2.67 x 10- 13 6^ 5 4 8 \ 5 «> 6 - 8 ^, if ^ = 0.32, ^ 
which indicates that the epoch when v a = v a > is very short, unless the ISM is very dense, 
e.g. n > 10 2 (Dai & Lu 1999, 2000). For the stellar wind case, the transition from v a = v a> 
to u a = v a< takes place when v a = v c , 



a 14 -0.80 0.07 f-l.02zr-0.44 4 1.02* :f f _fti 

aC ~~ 1 046f ~ a85 P - 08 ^-1-09/7-0.47 4 1.09/. 1f _ n 09 [ - b > 



The flux density at the observed frequency z/ from the synchrotron component for t < t ac 

is 

" (f) 2 (^)^, max oct( 1+fc+m )/( m+1 ), i/<i/ c , 

(-) 5/2 ^)- 1/2 ^ max OC t (8+*+^)/[4(m + l)] > ^ < „ < ^ 

(^-)- 1/2 F, max OC t (8-3*-4™)/[4 (m+ l)] > ^<^ m , 1 ^ 

(— r p/2 (^r i/2 i^ max oc t P-a-p(^» + fc)]/[4(m + i)] > ^ < V) 
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while for t ac < t < t cm , the flux density is 

l (|) 1/3 ^, max OC tdl-^Mm+l)] j U a <U<U e , 
( — )- p/2 (^)- 1/2 ^, max OC t [8-2*-p(4 ra+ *)]/[4(m + l)] > ^ < „ 



2.2. Constraint on the IC component in an X-ray afterglow 

The synchrotron-self-Compton (SSC) spectrum is featured by the characteristic IC fre- 
quencies, i.e., z/* c ~ 2^u a , v l £ ss 2^ m - m u m and z/* c ks 2^ 2 c v c . The IC frequency z/* c can be 
directly determined by 

vf = minH^,^}, (31) 

where v^ K « 2^ 1 c v a ^ < and i/^ ps 2^\v a> . Inserting equations (8), (16), (25) and (26) into 
the above equation, we obtain 



2.24 x 10»(1 + ^^Vs^^ Hz, ISM, 

3.01 x 10^1 + z)-XX°^f^ Hz, wind, 

(32) 

while the expression for v 1 ^ is 



,ic 



1.19 x 10^(1 + ,)-^ 2 f^!^^ Hz, ISM, 



>l/6 J ^cm,53"' V. 

4.27 x 10"(1 + ,)- 1 e;^L^ 1 2 k 1 % 6 < 5 3^*- 4/3 (7 L ) (5 - 14£)/[3(2 - e)1 Hz, wind. 

(33) 

As we can see, z/* c is below the X-ray frequency v ~ 10 18 Hz for typical parameters in most 
times during the fast-cooling phase. For simplicity, we do not consider this frequency for our 
estimation of the IC component in the X-ray light curve. 

The SSC frequency v 1 ^ equals to uf 1 when t = t cm , i.e. = 2^ cm u cm . According to 
equations (15) and (17), we obtain 

IC _ 1 , 3 5/4 m e C 3 / 2?77p >( 5fe _ 9 )/2 (l5fe-22)/8 (5fc-14)/8^(5fc-2)/4 

128 1 2tt J e 2 (l + z) l 3m e j e 

X ( i 443-^)-5/2 [ (3-fc)^ m ](5fc -6)/4 (34) 
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which can be further reduced numerically as 

= f !- 72 >< i 02 ^ 1 + ^XX'ZUyT E^n-^ Hz, ISM, 
""" ' 1.17xl0 19 (l + z)- 1 e ei _o.5e^ 2 2 . 5 Ci 2 /6^c m) 53^ 5/2 Hz, wind. 



The characteristic SSC frequencies z/* c and z/^ evolve with time as 

= J_)(2m+7fc-8)/[2(m+l)]^ ^IC = ^IC ( J_)-(6m+k)/[2(m+l)] _ (gg) 

The peak flux density of the SSC spectrum, -F*max> i s roughly the product of the peak 
flux density of the synchrotron spectrum by the Thomson optical depth. Considering some 
numerical factors of order unity, the exact expression is (Sari & Esin 2001) 

FlU - ^x (a T Rn)F„, max oc ^-^/p^i)] > (37) 

where xo ~ 0.5. The inverse Compton spectrum above z/* c is similar to the synchrotron one, 
which can be approximated by several power law segments, i.e. 

/ v Nl/3pIC .(16-llfc-4m)/[3(m+l)] IC 



" c 



,Tc 

c 

^_,-p/2 ( K^\-l/2 F lC .[8-2fe+4m-p(6m+fc)]/[4(m+l)] IC . 

jn / V . ,tc / 1 i/,max ^ ^ 5 ^ v i 



_^N-1/2 F IC f (8-3fc-2m)/[4(m+l)] IC < < IC 



where we have neglected the logarithmic term for v > v l £. 

The SSC flux density begins to dominate over that of the synchrotron radiation in the 
overall synchrotron + SSC spectrum at the crossing point, which corresponds to v~ (Sari 
& Esin 2001). Using equation (38) and the standard synchrotron spectrum, and assuming 
^x > v m > v c we obtain the crossing point frequency for two cases, i/* C < ^1° and z/* < 

, ,IC ^ ; .IC : „ 



, ^ C <--f[c 1 ^(^) 3p - 2 (2 7c 7e, min ) 2 - p ] 3/(2+3p) , if 4 C <^ C , 
' - I x C >-^[ci^(^) 3 - 2 (27 c 7e, nun ) 2 -1 1/( - 1) , if vf<vf<^, 1 J 

225(1 - e) 2 (p - l) 2 

where the coefficient is ci = tt. ; -777 777. To derive the above equation we have 

49xq(4 — k — e) 2 (p - 2) 2 4 

used the relation 

77.= 3(4-fc-6)(p-2) 

7c7e ' mm 8(1 - e)(p - 1)(1 + y) e B <7 T ni2" 1 j 
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Since l/(p — 1) is always larger than 3/(2 + 3p), one can determine v l £ directly by 

^ C = max{^ c < ,z/ I x c > }, (41) 
without judging whether v 1 ^ < z/* or not. 

We have calculated numerically the temporal evolution of v^p < and v l ^ > for typical 



physical parameters. In the ISM case, the expression for v 1 ^ < is 



( 3.5 x 10»(1 + z)-X^)-%E^n^^ Hz, p = 2.2, 

I 9.1 x 10"(1 + z)-X^t-L E ™>-™^ Hz > V = 2.4, 

(42) 

while the expression for i/^> is 

3.9 x io«(i + z )-^iS50- us (r)' &B/(HHz . ^ = 2 - 2 ' 

ic _ J V (43) 

3.0 x 1(F(1 + ^-^-i^c^^n-"^^)- 57 /^ Hz, p = 2.4. 



The crossing point frequency z/^ decreases throughout the fast cooling phase in the ISM case. 
The transition from the initial i/* c = ^x°> to the late z/^p = happens when u^p = z/^ c , 
i.e. at 

f fl O4-T-0-39 0.48 p0.04 „0.08^ ;f f _ n 1 

I U - ZUe e -0 .5 e B,-2.5- C/ cm,53 n f ™, II 6 — V.6Z, 

for p = 2.2, and 

IC _ f 0.057e-°_f e°A. 5 ^V°- 16 *«n, if c = 0.1, ( . 

\ 0.064e- a28 5 e0 J 4 l 2 . 5J B c m 06 53 n ai3 t cm , if e = 0.32, ^ 

for p = 2.4. The condition for the appearance of the IC component in the soft X-ray afterglow 
is that v^P a t tcm must be much less than v = 10 18 ui 8 Hz, which leads to the lower limit on 
the ambient density n. Using equation (42), we obtain the lower limit of n as 

f A Qfl _L -A-0.41.,-0.41 -1.27 -0.56 p-0.60 -3 _ 9 9 

~ 1 9^1+ -N-0.42 -0.42 -1.28 -0.56 p-0.60 -3 _ 9 4 

^ 2.0(1+2:; Z/ lg e e ,-0.5 e B,-2.5- C/ cm,53 Cm > P — 2.4. 

The lower limit of n for the emergence of IC component in the X-ray afterglow in the fast 
cooling phase is typically ~ 1-10 cm~ 3 (Panaitescu & Kumar 2000). 



In the stellar wind case, the expression for i/^< is 



4.5 x 10 18 (1 + ^)- 1 e^ 6 . 5 e^ 5 ^< 2 - 43 (^)( 3 - 1 - 5 - 7 ^[ 4 - 3 ( 2 -^ Hz, p = 2.2, 



1.4 x 10 18 (1 + ^-^^o.^^^^^ 2 - 3 ^^)^-^/! 8 ^ 2 -)] Hz, p = 2.4, 



(47) 
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while the expression for u^ > is 

IC j ^xlO^l + .r^e^ P = 2.2, 

1.1 x 10 17 (1 + ^)- 1 ee :T 5 ^ 2 .5^ m %^ 2 - 21 (T^) (7e - 26)/[7(2 - e)1 Hz, p = 2.4. 

6 r 



(48) 

The crossing point frequency decreases with i/* c = i/^> initially. However, the temporal 
behavior of the crossing point frequency at late times, v^P = ^x°o depends on e. The 
transition time when i/^ < = = z/ c c , is 



f n 7 r -0.26 0.22 zti-0.02 4 0.04/. jf ,_f)l 

t IC = J U -' de e,-l e B,-2.5- C/ cm,53 /1 * 'on, H t — U.I, 

I 0.56e e _' . 5 e^ 2 5 i? cmj53 A°- 04 t cm , if e = 0.32, 



for p = 2.2, and 



C n qq £ -0.28 £ 0.21 771-0.04 jo. 08/ ;f ,_ni 

.IC _ I U.JSCg^! e B -2.5-^0771,53 ^* r cm, 116 — U.i, . . 

x ' c ~ i n 97.-0-28 0.21 p-0.0440.08/ if , _ n 09 

I u - Z ' e e,-0.5 e B,-2.5- C/ cm,53 /1 * II 6 — U.OZ, 

for p = 2.4. After this time, the crossing point frequency z/Jp = v l £< wm increase for 
e < 31/57 (2/9), or continue to decrease for e > 31/57 (2/9) for p = 2.2 (p = 2.4). The 
emergence of the IC component in the soft X-ray afterglow requires that the minimum of 
z/^ during the fast cooling phase is below v = 10 18 z/i 8 Hz, which leads to a lower limit of A* 

as 

f 1.88(1 + .)- - 4 V-°- 41 6^ . 5 6^ 2 2 5 ^ 5 3, if 6 < 31/57, 

~ I 2.18(1 + z)^*U^™e£* B E%", jf £ > 31/57j ^ 



for p = 2.2, and 



r 0.94(1 + ^- - 43 <- 4 ^^V^-2.5«°53, if c < 2/9, 
' : ^' 1.15(1 + .)- - 42 <- 42 6°;\ 5 6^ 2 5 5 ^ 9 53 , if e > 2/9, 



(52) 



for p = 2.4. The above constraint on A* is insensitive to the other physical parameters. This 
implies that the contribution of the IC component to the X-ray afterglow can be neglected 
during the fast-cooling phase for typical values of A* < 1 as indicated from the observations 
of Wolf-Rayet stars and the fittings to some GRB afterglows (Chevalier & Li 1999; Chevalier, 
Li, & Fransson 2004). 



3. The late slow-cooling phase 



The hydrodynamic evolution of the blast wave in the slow cooling phase can be approx- 
imated as that in the early fast cooling phase. This approximation is validated by the fact 
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that the radiation efficiency of the blast wave in the slow cooling phase evolves as, 



6 = e e { — ) (p - 2)/ \ (53) 



v. 



which decreases very slowly with time as long as the electron power law index p does not 
deviate far from 2, e.g. p ~ 2.2 as expected in the relativistic shock acceleration theory 
(Achterberg et al. 2001). This fact will prolong the semi-radiative phase by at least two 
orders of magnitude in time than t cm for typical e e ~ 1/3 and p = 2.2 1 . Hereafter we assume 
e m e e for t > t cm . The time when e C e e happens is expected to be very late, when the 
Lorentz factor of GRB conical ejecta has already dropped below the inverse of the initial 
opening angle and the resulting light curves deviate from the spherical-like ones, which is 
beyond the scope of this paper. 



3.1. Properties of the synchrotron radiation 

The temporal behaviors of the typical frequency v m and peak flux density F v masL in the 
slow cooling phase are the same as in the fast cooling phase. However, the Compton param- 
eter in the slow cooling case is no longer a constant but evolves as Y = /_1(_!ZI)(p~ 2 )/ 4 . 

V € b v c 

Since the cooling Lorentz factor and the cooling frequency behave as 7 C oc (1 + Y)~ l and 
u c oc 7;? , we obtain 

_ /J_\(mp+4fc-4)/[2(4-p)(m+l)] 
Jc — )e,cm\, ) , 

v = v / J_\[6fc-8+(p-2)(4m+fc)]/[2(4-p)(m+l)] 
c cm \ , j 5 

Y = r±(_LyKp-2)(,rn+k-i)]/[(*-p)(™+i)]_ (54) 



e B t, 

The SSA frequency in the slow cooling phase is 



r c (p-l) enR l3/5 
V ^ Vm[ (3 - k) Bl lJ ' lf V * <V - 



Va ~ \ r cp(p-l) enR - - - (55) 

'Sa,> - Vm\ {3 _ k) B ^ 



/e.min 



^^The hydrodynamics will deviate significantly from that of the constant radiation efficiency e — e e in the 
very late times of the slow cooling phase when e < e~ x e e (e « 2.71828). This happens when t > 1300i cm 
(90t cm ) for the ISM (wind) case for p = 2.2 and e e = 0.32 (see equations 22 and 54). It will happen even 
later if we adopt the adiabatic relations for v c and v m . 
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which can be determined by v a = min{z/ a< , ^ a ,>} without judging whether u a < u m or not. 
The case for u a > u c can be neglected. The numerical expression for u a < is 

3.72«f (1 + ^-^-io^^^Cr/a^sa^^T^)- 367 ^ 4 -^ 1 GHz, ISM, 

16.4«f (1 + zre;XX%QT E2%m^r- m ^ )] GHz, wind, 

(56) 

where k p = (p — l){p + 2)/(p + 2/3), while the expression for v a>> is 

3.13 x 10"(1 + ^)- 1 e-i-§ 9 5 e^ 5 ^;f|n-°- 37 (-^)-( 8 - 6+e )/[ 3 - 1 ( 4 - e )] Hz, ISM, 
% 3.46 x 10"(1 + ^-^-i^e^L^^ - 3 ^^)^ 6 - 3 - 3 - 1 ^! 3 ^ 2 -^ Hz, wind, 



'on 



(57) 

for p = 2.2, and 

2.20 x 10"(1 + z)-^ 5 e^ 5 E^n^ M (^)-^y^^ Hz, ISM, 
] 2.89 x 10"(1 + z)-%^e^El^A^\^)-^^^ Hz, wind, 

tern 

(58) 

for p = 2.4. In the ISM case, the transition from the earlier v a — v a < to the later v a — v a > 
happens when v a — v m , at 

f IIyIOV'W ' 79 p-0.79 -1.38 / : f ,_ ni 

am ~ 1 Q ZL V in2 ff a38 f - a95 r - 76 p-0.76 T7 -1.33 / : f _ n on 

which indicates that the time when u a = i/ 0j> is very late, unless the medium is very dense. 
In the stellar wind case, the transition from u a = u a < to u a = u a > takes place when v a = u m , 
at 

f no 7 ^0.63 1.74 -0.16 /-2.22P0.95 4-2.22, - f , _ n. 1 

_ i vo.thup t e _ l t B _ 2 ^ 1 i & rj cmb?> i\ if L cm , 11 e — u.i, , . 

lam - S „„„ , 0.56 1.54 -0.14 ^1.96^0.84 ,4-1.96, : f _ n oo l DU J 

^ OOU.Ofi, p fc e ,-0.5 fc B,-2.5Sl/6 -^cra^^ 1 * L cm, u t — U.OZ,, 

which also indicates that the time when v a = is aj> is very late for the stellar wind case. 
Therefore we neglect the case of t > t am below. 

The flux density at the observed frequency v from the synchrotron component for t cm < 

) ^ ^i/,max (X t ^ \ V < Z/ a , 



(^) 1/3 ^, m ax OC t (»-«-™)/P(«H-D] , „ a < „ < „„,, 

( — )- (p - 1)/2 ^, m ax OC t (12-5,-p(4 m+fe ))/[4( m+ l)] ; I/ m < 1/ < I/„ 

(_)-P/ 2 (^)-(p- 1 )/ 2 i? ymax (X i {(4-p)[12-6fc-4m-p(4m+fc)]+8(fc+m-l)}/[4(4-p)(m+l)] ^ < y 

(61) 
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3.2. Constraint on the IC component in an X-ray afterglow 

The temporal behaviors of the typical SSC frequency v£ ~ 2 r y1 niin v m and peak flux 
density i^^ax m the IC spectrum are the same in the slow cooling phase as in the fast 
cooling phase. The IC frequency corresponding to v c is 

V \Q ^ 2 ^jy c = z/ IC^^^(6mp+pfe+12fe-8m-16)/[2(4-p)(m+l)]^ ^ 



The IC frequency z/* c can be directly determined by z/* c = minjz/^, z/^}, where ~ 
2 7e,min i/ a,< and 1/*°. w 2^ min i/ aj> . Inserting equations (6), (16), and (56) - (58) into the 
above equation, we obtain 



1.75x lO^l^)- 1 ^^ Hz, ISM, 

2.36 x lOV 5 ^*)- 1 ^'^^ Hz, wind, 



V a,> = 



(63) 

while the expression for v a> is 

IC _ } 1-48 x 10»(1 + .)- 1 6-io 94 5^ 2 °5^;I^ ISM, 

4.99 xlO"^*)" 1 ^^^ Hz, wind, 

(64) 

for p = 2.2, and 

IC j l^xlO^l + .r^e^ ISM, 
9.35 xlO 17 ^,)- 1 ^.^^ Hz, wind, 



"a,> 



"cm 



(65) 

for p = 2.4. As we can see, z/* c is below the X-ray frequency v ~ 10 18 Hz for typical 
parameters. We thus do not consider this frequency in our estimation of the IC component 
in the X-ray light curve in the slow-cooling phase. 

The inverse Compton spectrum is 

'— W 3 F IC oc f J_^(12-7fc)/[3(m+l)] ic 
„IC' v,max ^ V . ) i u ^ u m i 

v ra l cm 

■ _^_\-{p-l)/2 F IC ( J_ ) [W-9k+2m-p(6m+k)]/[4(m+l)] IC IC 

'J^\-p/2/^c_\-(p-l)/2rfC / * \[6(2-fc)(4-p)+p(6mp+pfc-20m-4)]/[4(4-p)(m+l)] IC ^ 

(66) 



,iIC 
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where we have neglected the logarithmic term for v > z/* c and also do not consider the lowest 
spectral segment below z/* c for simplicity. The relation between the peak flux density of the 
SSC spectral component and that of the synchrotron component is (Sari & Esin 2001) 

(p-l)(p+l/3) . T-J \ T-I 

The critical frequency corresponding to the crossing point of the synchrotron spectral com- 
ponent and the SSC component is 

^x C <-^ C [c 2 -( 7 ^) 4 (27 C 7e, ml n) 2 -1 3/(2+3p) , if ^ X C <" IC 



ic _ J e e 7 e 



..IC - ,,IC[„ 6 £/ 7c \4/r. \2-pl :f ..IC ^ ..IC ^ ..IC 1 ' 

^x,> = [ C 2 — I" J (/7c7e,minj J, U ^ m < ^ x < > 

7e,min 



where we indude the coefficient e 2 = -^pt^M^L-, which is much larger 

93Jq(4 — k — c) (p — 2) (p + 1/3) 

than unity (by at least one order of magnitude), but was neglected in equation (5.1) of Sari 

3 

& Esin (2001). Since is always smaller than unity, i/* c can be determined directly by 

2 -|- 3p 

v 1 ^ = max{i/Jp ^x C >} without judging whether i/* c < or not. 

We have numerically calculated the temporal evolution of and ^x°>- ^ n ^ e 
case, the expression for v!P < is 



IC j 3.5 x 10»(1 + s)"^^^^^ Hz,p = 2.2, 

7.3 x 10 18 (1 + zr^^e^E^n-^S^-^y^^ Hz, p = 2.4, 



1 cm 



(69) 



while the expression for v l £ > is 



IC j 1-7 x 10"(1 + ,)- 1 6-% 7 5 e^, 5 ^>- 2 ^ p = 2.2, 

1.8 x 10"(1 + ^"^^o" c^S^- 2 ^^)^^ 4 -) Hz, P = 2.4. 



(70) 

We can see that v 1 ^ decreases first with v~ = <, then increases with i/* c = >■ The 
time when u 1 ^ reaches its minimum, v l ^ < = v\ p> = is 

C zL -0.47 p-0.04 -0.08/ ;f f _ n 1 

/IC ] ^• 0fc e,-l fc B,-2.5- C/ cm,53' i 'cmi lie — U.l, . . 

x ' m I c; 1.0.34 -0.41 p-0.04 -0.07/ :c , _ n qo I ' ^ 

I d - ifc e,-0.5 fc _B,-2.5- C/ cm,53' i 'cmi He — U.OZ, 



for p = 2.2, and 



f 1 7 R c 0.33 ,-0.48 1^-0.07^-0.15/ ;f , _ n 1 

/IC _ J 1 ' - De e ,-l e B,-2.5- C/ cm,53 n ^cm, 116 — U.l, 

l x,m — 1 1« 7^0.29 ,-0.41 p-0.06_-0.13/ :f , _ A 19 ^ ' 

I ib - 7e e -0.5 e S ,-2.5 jfc/ cm,53 n *cro, « C - 0.32, 
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for p = 2.4. The IC component could appear in the X-ray afterglow only if the minimum of 
z/^ c is less than v = 10 18 ^ 18 Hz, which leads to a lower limit on the ambient density n. We 
obtain the lower limit of n as 



> 



for p = 2.2, and 



> 



8.6(1 + z^u^e^e^E^ cm"', if e = 0.1, 
1.6(1 + z^v^e-l&e^EZ™ cm" 3 , if e = 0.32, 



1.9(1 + 2)^%^^ cm- 3 , if e = 0.1, 
0.5(1 + z)^«v^el™e^E£$ cm" 3 , if e = 0.32, 



(73) 



(74) 



for p = 2.4. This lower limit of n for the emergence of IC component in the X-ray afterglow in 
the slow cooling phase is typically in the range of 1 — 10 cm~ 3 (Sari & Esin 2001; Panaitescu 
& Kumar 2000; Zhang & Meszaros 2001). However, the true lower limit of n is even smaller 
than that given in the above equations, since we have neglected the case of u 1 ^ > z/ c '. The 
spectral segment when i/Jp > v 1 ^ is oversimplified by a single power law approximation. In 
fact, the logarithmic term dominates at higher frequencies. The true evolution of v 1 ^ is 
always decreasing, although the decreasing rate is slowed at late times (Sari & Esin 2001). 

In the stellar wind case, the expression for i/^< is 

IC _ } ^xlO^^,)-^;^ p=2.2, 
3.5 x 10 18 (1 + ^^^^o^^^^^^^X)-^^)]/^^-,)] Hz p = 2A 

(75) 

while the expression for i/^> is 

IC _ } 7.2 xlO"^*)"^^^ P=2.2, 
3.0 x 10 16 (1 + .)- 1 e -°jL^ 2 .5«,53< 2 - 1 (7^) (5 - 4 " 3 - 5e)/(2 - e) Hz, p = 2.4. 

"cm. 



The time when z/^ c reaches its minimum, i/^< = i/^> — v m-> ls 



(76) 



14,0.29,-0.24 po.02 4-0.05^ :f , _ n 1 

JC ^_ J i -^ e e,-l e B,-2.5- C/ cm,53 /1 * l cm, II £ — U.l, . . 

L x,m 1 I q,0.29 c -0.24 77.0.02 a -0.05+ ;f , _ n 00 ^ ' 

1 - yt e,-0.5 t _B,-2.5- C/ cm,53^ 1 * ^m, 11 t — U.OZ,, 



for p = 2.2, and 



q 1 ,0.32 -0.23 p0.04 a -0.09+ :f , _ n 1 

-IC J " 3 - ie e,-l e B,-2.5 jC/ cm,53 /i * r cm, II t — U.l, . . 

r x- m 1 44.O.32 -0.23 £0.04 A-0.08+ if f - 32 1 } 

^•^ fc e,-0.5 fc _B,-2.5- C/ cm,53 /1 * L cm, 11 t — U.OZ, 
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for p = 2.4. The emergence of the IC component in the X-ray afterglow requires that the 
minimum of u 1 ^ is lower than the X-ray frequency v = 10 18 z/ 18 Hz, which leads to a lower 
limit on A*, i.e. 



1.32(1 + zY^r^E^ if e = 0.1, 



~ 1 1.55(1 + .)-°- 4 Vr 8 - 42 e^ 5 . 5e ^ 2 . 5 ^ 53 , if e = 0.32, 



(79) 



for p = 2.2, and 



> 



1.03(1 + :) iU:i /^ 1 - ;:1 l : i: ;/^;;;l : , if e = 0.1, 
1.13(1 + ^)" - 43 ^r8 a43 <^o.5« 9 53, if e = 0.32, 



(80) 



for p = 2.4. The above constraint on A* is insensitive to other physical parameters. Together 
with the constraint on A* for t < t cm , we conclude that the contribution of the IC component 
to the X-ray afterglow is insignificant and can be neglected for A* < 1 as indicated from 
observations of Wolf-Rayet stars and fittings to some GRB afterglows (Chevalier, Li, & 
Fransson 2004; Panaitescu & Kumar 2001, 2002). 



4. Afterglow light curves of semi-radiative blast waves 

We assume below that the physical parameters do not deviate significantly from that 
chosen in previous sections. The contamination of the IC component in the high frequency 
afterglow (e.g. the soft X-ray afterglow) is not considered for simplicity. However, the IC 
emissions can be inferred from equations (38) and (66). Under these assumptions, the light 
curve at an observing frequency v can be determined by comparing the frequency with the 
critical frequencies v cm and u ac , where v ac is the SSA/cooling frequency at t ac and can be 
calculated from equations (22), (27) and (28). Roughly, there are three types of afterglow 
light curves in various frequency ranges separated by these two critical frequencies. A careful 
inspection of the order of the transition time t cm and the crossing times t a , t c , t m gives four 
types of light curves for both the ISM case and the stellar wind case. The crossing times t a , 
t c and t m correspond to the times when the frequencies u a , u c and u m equals the observing 
frequency, respectively. 

In the case of ISM, the orders of these times are (A) t c < t a < t m < t cm for v > u ac ; 
(B) t a < t c < t m < t cm for v cm < v < u ac ; (C) t a < t cm < t m < t c for v a (t cm ) < v < u cm ; 
(D) t cm < t a < t m < t c for v < v a (t cm ). Here v a (t cm ) is the SSA frequency at t cm . Using 
the equations (29), (30), and (61) in the case of k = 0, the light curves in each case can be 
constructed as (A) F v oc tV (t < t c ), F v oc 4(5-20/(4-0^/2 ^ f ^ k t - (i+2 e )/(4- e ) u -i/2 

(t a , t m ), F v OC r (3p-2+20/(4-0 i/ - P /2 (tm; tcm); Fu r (3 P -2 + 2e)/(4-e) + (2 + e)(p-2)/(4 + e)/(4- p)i/ - P /2 
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(t > t cm ); (B) F v oc fV (f < t a ), F v oc ^-1^/3(4-^1/3 F|/ oc r (1^/(4-^^-1/2 

(f c , t m ), F„ OC t"(3p-2+2 e )/(4- e ) z/ -p/2 j \ ^ K t _(3p_2+2e)/(4- e )+(2+ e )(p-2)/(4+ e )/(4-p) I/ -p/2 

(f > L!)J (C) ^ OC fV (f < t a ), F„ k'^H.I/S ( to> tcm ), J?„ oc ^2-3^/(4-^)^1/3 ( ^ 
t m ), F v OC t-3(p-l+^)/(4-e) z/ -(p-l)/2 U p v K i -(3p-2+2e)/(4-e)+(2+e)(p-2)/(4+e)/(4-p) l/ -p/2 

(t > t c ); (D) F„ oc tV (t < t cm ), F v Kt^ 1 ^-^ (t cm , t a ), F v oc f (2-3 e )/(4- e)l/ i/3 ^ 

F„ OC t-3(p-l+ £ )/(4-,) z/ -( P -l)/2 (tm; ^ ^ r (3 P -2 + 2e)/(4-e)+(2+e)(p-2)/(4+e)/(4- p)i/ -p/2 ^ > g 

The light curves in the ISM case are illustrated in Figure 1. The crossing times t c and t a 
in case A and t a in case B occur very early in high observing frequencies, while t c in case 
C and t a , t m and t c in case D occur very late in low observing frequencies. We thus neglect 
these crossing times in the figure. As indicated in Figure 1, the radiation efficiency, e, has 
a marked effect on the afterglow light curves. It changes the temporal decaying index a 
(defined as F v oc t~ a ) of the light curve significantly. For illustration, we adopt e = e e = 1/3 
and p = 2.2 in the following. The initial slowly increasing light curve segment, F v oc t 1 / 6 , 
predicted in the standard adiabatic blast wave model will changes to be a slowly decreasing 
one, F v oc t( 2 ~ lle )/ 3 ( 4 - £ ) t -0 - 15 , as shown in the early segment of case C. This makes the sub- 
millimeter afterglow less competitive to distinguish between the ISM and the stellar wind, 
as proposed by Panaitescu & Kumar (2000), if the observations are not frequent enough. At 
the optical wavelength, the early light curve behaves typically as F v oc t _ ( 1+2e )/( 4 - e ) ~ t -0 - 45 
rather than F v oc t -1 / 4 in the adiabatic case. When v m crosses the observing frequency, 
i.e. t > t m , the optical and X-ray light curves decays as F v oc £-(3p-2+2e)/(4-e) ^ t~ L44 , as 
shown in cases B and A. It should be noted that many X-ray afterglow light curves and a 
considerable fraction of optical afterglow light curves have temporal decaying indices steeper 
than predicted (3p — 2)/4 ~ 1.15 by the standard adiabatic model (for p ps 2.2). The 
observed decaying indices in the X-ray afterglow light curves are (ax) = 1-33 ± 0.38, while 
the median value of the observed X-ray spectral indices f3x, F v oc z/ _/3x , is ~ 1.05 (Berger, 
Kulkarni, & Frail 2003; De Pasquale et al. 2003). Assuming the X-ray frequency is above 
the cooling frequency and the spectral index is f3x — p/2, the measured p is consistent with 
the standard value of the index of electron energy distribution, p = 2.2 — 2.3, predicted in the 
relativistic shock acceleration mechanism (see Achterberg et al. 2001 and reference therein). 
However, the observed mean temporal decaying index (ax) requires a relatively larger (p), 
~ 2.44, provided the shock is adiabatic. There are several caveats on the observations of the 
X-ray afterglows. First, the temporal behavior of the X-ray afterglow is hardly influenced 
by the equal arrival time surface effect, which will mix the earlier light from high latitudes 
into the present light. This effect is important especially for high observing frequencies, e.g. 
the optical and X-ray. The profile of surface emissivity of the relativistic shock is ring-like in 
these high frequencies (Sari 1998). This will moderately slow down the decreasing rate of the 
afterglow after t m for the theoretical light curve A. However, the X-ray afterglow is immune 
to this effect because its t m is very early and the observed decreasing index ax is based on 
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the observations typically several hours after the burst. Second, the measured fix is reliable 
since the X-ray absorption in the medium along the line of sight takes place at v < 1 keV 
while the observing window is ~ 2 — 10 keV. Lastly, one should be cautious when interpreting 
the property of X-ray afterglow with the synchrotron radiation mechanism, because the X- 
ray afterglow may be contaminated by the synchrotron-self- Compton component. However, 
there have so far been only a few X-ray afterglows that were confirmed to have the IC 
components. Therefore, the radiative corrections to the afterglow light curves must be taken 
into account based on the observations. As can be seen in Figure 1, the light curve at 
high frequency (type A and B) flattens when the afterglow enters the slow-cooling phase, 
t > t cm . At the transition time t cm , the spectrum nearby the observing frequency changes 
from v c < v m < v to v m < v c < v, while the expressions for the flux density are the same, 
i.e. F v = -F^max^m v x J 2 v^ v I 2 . The flattening of the light curve results from the Compton 
parameter Y in the flux density, i.e. F v oc Y -1 , since v c oc (1 + Y)~ 2 m Y~ 2 . The Compton 
parameter Y in the slow-cooling phase decreases slowly, contrary to its constancy in the 
earlier fast-cooling phase. From equation (54) for Y in the slow-cooling phase and adopting 
k — 0, the change of the temporal index around t cm is Aa = (2 + e)(p — 2)/[(4 — e)(4 — p)\, 
which is shown in the last segments of panels A and B in Figure 1. Note that since Sari et 
al. (1998) did not discuss IC cooling, there is no related segment in their e = light curves. 

In the case of stellar wind, the orders of the crossing times are (A) t a < t m < t cm < t c for 
v > v cm ; (B) t a <t c < t cm < t m for v ac < v < v cm ; (C) t c <t a < t cm < t m for v a (t cm ) < v < 
u ac ; (D) t c < t cm < t a < t m for v < u a {t cm ). Using the equations (29), (30), and (61) in the 
case of k = 2, the light curves in each case can be constructed as (A) F v oc £( 7 - 5 0/ 2 ( 2 -0 2,5/2 
it < t a ), F v oc t-WP-Ou-V* (t a , t m ), F u oc (t m , t cm ), F v oc 

r [3 P -2-(p-2) e ]/2(2- e )+(p-2)/(4- P ) I/ - P /2 Fy K f - [3p-l-(p-l) e ]/2(2- e ) I/ -(p-l)/2 ^ > Q. (g) 

F„ OC ^7-50/2(2-0^/2 (t < Fy K r (l +e) /2(2-e) i/ -l/2 ^ ^ K r (4-^/3(2-^1/3 

t cm ), F v OC r^^l/VS (tcmj tm)j K K r[ 3p-l-(p-l)e]/2(2-e) i/ -( P -l)/2 (f > ^ (C) ^ K 
^7-50/2(2-0^5/2 (t < tc)j ^ K ^4-30/(2-0^2 ^ ^ ^ r (4- e )/S(2- e)l/ l/S ^ > ^ K 

r 03(2- e)i/ l/3 (fcm> tm); ^ K r[ 3p-l-(p-l)e]/2(2-e) i/ -(p-l)/2 (f > ^ . (D) ^ K ^7-50/2(2-0^5/2 
(f < t e ), F v OC ^(4-30/(2-0^2 (tc; Fv K t 2(l- e )/(2- e)l/ 2 ^ ^ K ^/Sp-^l/S 

f m )» oc t-[3p-i-(p-i)^]/2(2-0 z/ -(P-i)/2 ( t > tm y The light curves - m t he stellar wind case 
are illustrated in Figure 2. The crossing time t c in cases C and D occurs very early at 
low observing frequency, while t m in case D occurs very late. We neglect these crossing 
times in this figure. The radiation efficiency has a significant effect on the light curves in 
the wind case. The flux density in the optical/infrared light curve decays initially with 
£-(!+0/ 2 (2-e) ^ £-0-4^ raTj h er than t -1 / 4 in the adiabatic case, which is shown in case A. For 
the X-ray afterglow the crossing time t m when the typical frequency u m crosses the observing 
frequency is much earlier, the light curve behaves as i _ [ 3 P _2 -(p-2)e]/2(2-e) ^ £-i-36 during the 
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whole fast-cooling phase, which is more consistent with observations than the adiabatic light 
curve (Berger, Kulkarni, & Frail 2003). The light curve at high frequency flattens when the 
afterglow transits to the slow-cooling phase. By the same way as in the ISM case, from 
equation (54) for Y and adopting k = 2, the change of the temporal index around t cm is 
Aa = (p — 2)/(4 — p) in the wind case. Since Chevalier & Li (2000) did not include IC 
cooling, there is no relevant flattening segment in their e = light curves. Although the 
flattening of the optical/X-ray light curve around t cm predicted in the inverse Compton 
dominated cooling regime in the stellar wind case is more obvious than in the ISM case, the 
change of the temporal decaying index is only Aa ~ 0.1 for the former case. The detailed 
theoretical optical light curves taking into account the equal arrival time surface effect and 
the large error bars in X-ray afterglow observations prevent us from the identifications of 
such flattening. 



5. Conclusions and Discussion 

In this paper, we have investigated analytically the GRB afterglow hydrodynamics and 
constructed the semi-radiative light curves realistically. We focus on the case that the elec- 
tron cooling is in the inverse Compton dominated regime, i.e. e e > eg or 7 > 1. The 
realistic hydrodynamics is applicable for spherical blast waves with the assumption that 
the electron energy equipartition factor e e is not much larger than 1/3, which seems to be 
reasonable based on theoretical expectations and observations as well. In fact, the analyti- 
cal solution for afterglow hydrodynamics is almost tenable throughout the relativistic stage 
when e e < 2/3 (see equation 2). The only uncertainty is the actual evolution of the radiation 
efficiency in the late slow-cooling phase. Given p ~ 2.2 — 2.3, we conclude that a constant 
radiation efficiency is a good approximation for a fairly long time in the slow-cooling phase. 
The transition from fast-cooling to slow-cooling happens much later in the IC dominated 
cooling regime than in the purely synchrotron cooling regime. Since the actual radiation 
efficiency decreases very slowly in the slow-cooling phase, the semi-radiative epoch is further 
prolonged typically by two orders of magnitude in time, other than the whole fast-cooling 
phase as commonly used. As the GRB ejecta sweeps up more and more external medium, the 
Lorentz factor 7 of the shock decreases. When 7 equals the inverse of the initial half-opening 
angle 8 of the GRB conical ejecta, or jet, the following afterglow light curves deviate from 
the previously spherical-like ones, and the light curves in this paper are not suitable at this 
late stage. Generally, our semi-radiative afterglow light curves hold true when the shock is 

relativistic after the initial deceleration time and before the jet-like stage, with the radiation 

2 2 

efficiency satisfying e « e e <max{ — — -, -}. 



-22 - 



The adiabatic afterglow light curves in the synchrotron dominated cooling regime have 
been well studied in previous works (Sari et al. 1998; Chevalier & Li 2000; Granot & Sari 
2002). Sari et al. (1998) have considered the light curves of a fully radiative (e = 1) blast 
wave. However, our analytical results can not be directly applied to this case. Actually, in 
the fully radiative case, equation (2) can also be directly integrated and the hydrodynamics 
and the resulting light curves are the same as derived by Sari et al. (1998). The radiative 
corrections to the afterglow light curves in the wind model have been discussed by Chevalier 
& Li (2000). They had adopted the scaling laws of a semi-radiative shock given by Bottcher & 
Dermer (2000). The temporal exponents of the hydrodynamics and light curves in Bottcher 
& Dermer (2000) differ from ours. In fact, the radiation corrections for these exponents 
in their work are smaller than that in our work 2 . Cohen, Piran & Sari (1998) had studied 
the hydrodynamics of a semi-radiative relativistic blast wave considering a more complicated 
post-shock material distribution (similar to Blandford & McKee 1976) than that of the simple 
thin-shell approximation as in our work. However, the hydrodynamic evolution differs in all 
these three treatments of Cohen et al. (1998), Bottcher & Dermer (2000) and this paper. 
The hydrodynamic self-similarity index m we obtain lies between the other two's. Despite 
the difference between the hydrodynamics due to different approximations/treatments, we 
can see that the radiative corrections to afterglow light curves are significant. For example, 
the temporal decaying index of X-ray afterglows at around 10 hours since the main bursts 
varies in the range of 1.23 — 1.69 in the ISM case, and 1.21 — 1.58 in the stellar wind case 
for p = 2.2 — 2.3 and e e = 0.1 — 0.5, provided that the IC component is neglected. This 
range is consistent with the observations (Berger, Kulkarni & Frail 2003). In observations, 
the value of e e and p can be inferred for a particular X-ray afterglow from the so called 
"closure relation" between the temporal index ax and spectral index 0x, provided that 
the IC component can be neglected. Such a closure relation can be obtained from case 
A described in §4 in either ISM or wind cases. The application of this method to optical 
afterglows should be cautious. The early optical light curve has a broken power law profile 
around t m . In contrary to the early X-ray afterglow whose t m is much earlier and which can 
be regarded as a simple power law one, the optical light curve might be affected moderately 
when the equal arrival time surface effect is included. Another reason is that the reddening 
in optical spectra is unknown. While the Galactic extinction can be empirically decoupled, 
the extinction and reddening within the circum-burst environment and host galaxy are less 
constrained. These two facts prevent the credible measurements of the temporal and spectral 
indices of optical afterglows, respectively. 



2 The coefficient of the e term in these temporal exponents in their work is smaller than ours by an exact 
factor of 2. 
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We have got the criteria for the emergence of IC components in soft X-ray afterglows, 
through giving the lower limits of external medium densities. In the ISM case, the lower limit 
of n is ~ 10 cm -3 in the fast-cooling phase, while it is ~ 1 cm -3 in the slow-cooling phase 
(Panaitescu & Kumar 2000; Sari & Esin 2001; Zhang & Meszaros 2001). These are typical 
densities of interstellar media in our galaxy. In the wind case, the lower limit of the wind 
parameter is always A* ~ 1 (Panaitescu & Kumar 2000). Such a critical A* is also typical for 
Wolf-Rayet stars in our galaxy. It should be noted that the wind parameter obtained from 
fitting afterglows within the wind interaction model seems to be quite small (Chevalier & Li 
1999, 2000; Panaitescu & Kumar 2001; Dai & Wu 2003; Chevalier, Li, & Fransson 2004). The 
contradiction may be due to the limitation of our knowledge about the mass losses of massive 
stars at the last stage before their collapses. Taking the inferred low A* from afterglows, we 
draw a tentative conclusion that the IC components in X-ray afterglows are insignificant in 
the stellar wind case. For such a low A*, a wind bubble would be produced and surrounded 
by either an outer giant molecular cloud, a slow wind at previous evolutionary stage, or 
an extremely high pressure in a star-burst environment (Dai & Wu 2003; Chevalier, Li, & 
Fransson 2004). The termination shock radius of the wind bubble will be reached by the 
GRB shock within hours in observer's frame. The environment of the shock will change from 
wind type to uniform medium type after this time. We have neglected such a complicated 
case in this work. Recently, Yost et al. (2003) had relaxed the assumption of a constant 
magnetic equipartition factor e# and broadened the circum-burst medium types. They had 
made a detailed comparison between the results of different assumptions as well as different 
circum-burst media, and found the degeneracy of different assumed evolutions of eg and 
different medium types. Anyway, we adopt the constant e B assumption and consider the 
most possible medium types, i.e. the ISM and the stellar wind. 

The radiative corrections in modelling afterglows are important. Such an effect should 
be taken into account seriously in analysis of afterglows, especially when a large number of 
afterglows will be observed in the Swift era. It will affect directly the energetics of GRB 
remnants and therefore the actual efficiency of prompt GRBs (Wu et al. in preparation, 
which also includes the ICS effect). 
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Fig. 1. — Characteristic afterglow light curves from a semi-radiative blast wave in an in- 
terstellar medium (ISM) in various frequency ranges. The panels are ordered from high 
frequency (type A) to low frequency (type D). X-ray and optical light curves are typically of 
types A and B, while sub-millimeter and radio light curves are typically of types C and D. 
The physical parameters for these light curves are E cm = 10 53 ergs, n — 1 cm -3 , e = e e = 1/3, 
e_e = 10~ 2 - 5 , p = 2.2. The event is assumed to be at redshift of z — 1. 
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Fig. 2. — Characteristic afterglow light curves from a semi-radiative blast wave in a stellar 
wind in various frequency ranges. The panels are ordered from high frequency (type A) 
to low frequency (type D). X-ray, optical, and infrared light curves are typically of type 
A, sub-millimeter light curve is typically of type B, while radio light curve is typically of 
types C and D. The physical parameters for these light curves are E cm = 10 53 ergs, A* = 1, 
e = e e = 1/3, e B = 1CT 2 - 5 , p = 2.2. The event is assumed to be at redshift of z — 1. 



